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Abstract. Operationally, index functions of variable Hilbert scales can be 
viewed as generators for families of spaces and norms. Using a one parameter 
family of index functions based on the dilations of a given index function, a new 
class of scales (dilational Hilbert scales (DHS)) is derived which generates new 
interpolatory inequalities (dilational interpolatory inequalities (DII)) which 
have the ordinary Hilbert scales (OHS) interpolatory inequalities as special 
cases. They therefore represent a one-parameter family generalization of OHS, 
and are a precise and concise subset of VHS approriate for deriving error 
estimates for deconvolution. The role of the Hilbert scales in deriving error 
estimates for the approximate solution of inverse problems is discussed along 
with an application of DHS to deconvolution sharpening. 

1. Introduction 

In the analysis of the numerical performance of traditional regularization meth- 
ods (Engl et al. [1]), interpolatory inequalities between the norms \\u\\ a , a € R+, 
generated by an appropriate family of Hilbert spaces H a = R(T a / 2 ), play the cen- 
tral roleQ In terms of the original concept of a Hilbert scale, as introduced by Krein 
and Petunin [TO] and generated by a densely-defined, unbounded, self-adjoint and 
strictly positive operator T, such inequalities take, for a given value 7 of the linear 
scale between given values a and /3, the form 

(1) \\u\\ eaH i-e)p < HI* \\u\\y, for all #e (0,1) and < a < /3 

where (•, •) and |j • || denote the inner product and norm in the original Hilbert space 
H = Hq and where 

= ||T Q/2 u||, a>0 

and similar for 1 1 it || ,3. 

Through the appropriate choice of T and the values for a and b, the corresponding 
inequality (1) can be used to derive error estimates for the regularized solution of 
improperly posed operator equations that simultaneously take account of both the 
compact and smoothing nature of the operator (Groetsch [5], Natterer |17| . Schroter 
and Tautenhahn [15], Tautenhahn [2"U]). 

When utilizing such inequalities to derive error estimates for linear improperly 
posed problems, it was observed by various authors that realistic error estimates 
could only be derived for a subset of linear improperly posed operator equations. 
This led to the need to construct more general counterparts of the inequality (1). 
Hegland [Sj , [7j , by exploiting the spectral decomposition of an appropriately chosen 
operator T, first introduced the concept of a variable Hilbert scale (VHS) for a quite 
general index function <j). Then, by invoking appropriate regularity about the choice 
of the index function, Hegland established how more general counterparts of the 
interpolatory inequality (1) could be constructed. 



Here R(T°"/ 2 ) denotes the range of the operator T a / 2 . 
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The utility of this basic concept of a VHS, in deriving more representative error 
estimates and convergence rates for the regularized solution of improperly posed 
operator equations, has been subsequently exploited and/or modified by various 
authors including Nair et al. [16] , Mathe and Pereverez [12], [15] . Hofmann and 
Yamamoto [9], Mathe and Tautenhahn [T3] and Begout and Soria [2]. 

The concept of a variable Hilbert scale not only overcomes the mentioned short- 
coming of interpolatory inequalities of the form (1), but also allows a greater variety 
of interpolatory inequalities to be constructed. It is this aspect that is pursued in 
this paper. A recent example can be found in Nair et al. 16]. In our paper, the 
concept of a family of variable Hilbert scales and associated index functions are 
introduced using the spectral theorem for positive definite self-adjoint operators. 
The special subclass of dilational Hilbert scales (DHS) is then defined, and a new 
type of interpolatory inequality is derived and applied. 

The motivation is the analysis of deconvolution by sharpening (Hegland and 
Anderssen [5]). When the Gaussian model 

G 7 (f) = -G{t/i), G(t)=exp(-t 2 /2), 
7 

defines the smoothing kernel B 1 (e.g. the broadening function that models the 
measurement process in spectroscopy as examined in Section 3.3) in the convolution 
equation 

(2) 9 = B 1 * /, 
the convolutional relationship 

G a * Gp — G 7 , 
can be used to factor 

9 = G 1 * / 

in the following manner 

g = G a *Gp* f = G a * z, z = G *f 

with a determining the nature of the sharpening to be performed and (3 the form of 
the source condition. In this way, a and (3 perform a trade-off, mentioned earlier, 
between the sharpening and the achievable rates of convergence. Though such a 
strategy does not in general hold for non-Gaussian models of the smoothing kernels, 
it can be adapted to dilational models 

B~(t) — —B(t/j), Bit) ~ some representative peak function, 

7 

by performing the factorization in the following manner 

(3) g = B J:j3 * Bp* f = B lt p* z, B y = B yi0 *Bp, z = Bp*f, 

where the form taken by B ly p{t) depends on i? 7 , chosen to reflect the convolutional 
smoothing that the measurement system has performed on broadening / to give g, 
and on /3, chosen to ensure that the source condition encapsulated in z = Bp * / 
yields acceptable convergence rates. A discussion about how y — Bp * / provides 
an appropriate source condition is given in Section 13.31 

Now, however, it is necessary to determine and work with B lt p(t) as the sharp- 
ening to be utilized when performing deconvolution. The value of 7 reflects the 
difference in instruments with small 7 corresponding to expensive accurate mea- 
surement instrumentation. The scaling I/7, introduced into the definition of By, 
is an essential regularity condition. As (3 tends to zero, Bp(t) must have a be- 
haviour like that of the identity operator so that, under appropriate circumstances, 
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z recovers /. For the scaling chosen here, it guarantees that z — f when / is a 
constant. 

In recording a spectrum, two separate independent steps are involved: 

(i) sample preparation, and 

(ii) measurement of the spectral response of the sample. 

Here, equation ([2]) is a model for (ii). For (i), a separate model is required which, 
if convolutional, would take the form 

f = S*f 

where / denotes the ideal form of the sample and S the smoothing resulting from 
the sample preparation. 

In monodisperse molecular weight distribution modelling in rheology, / is a Dirac 
delta function of a single molecular weight [H] . In NIR spectroscopy, / is a collec- 
tion of particles of a single size and homogenously arranged [5] . In mass spectro- 
scopic peptide finger printing, / is a sum of Dirac delta function molecular weights 
corresponding to a perfect segmentation of the protein into individual molecular 
weight peptides. 

For the analysis in the sequel, this relationship between / and / justifies the 
conclusion below that / is an L2 function. If it was assumed that / was a Dirac 
delta function, the analysis below would have to be performed with weaker norms. 

The paper has been organized in the following manner. The concept of a DHS 
is introduced in Section 2 as a special case of a variable Hilbert scale (VHS). The 
generalization of Begout and Soria [5] of the Holder inequality for LP spaces is first 
used in Section 2 to derive a general form for variable Hilbert scales interpolatory 
inequalities, which is then utilized to derive the basic DHS interpolation inequality. 
The application of DHS in the construction for error estimates for the approximate 
solution of inverse operator equation problem is examined in Section 3 along with 
a discussion of analytic differentiation, source conditions and the application of the 
results to deconvolutional sharpening. 

2. Variable and dilational Hilbert scales 

The spectral theorem for a positive definite selfadjoint operator T on a Hilbert 
space H (see, e.g., [H]) introduces (for each T) a family E(X) of orthogonal pro- 
jections such that Ef :9 (X) := (E(X)f,g) defines a Stieltjes measure on (0,oo) and 

/>oo 

(Tf,g)= / XdE f JX) 
Jq 

for all g € H and / S D(T) (the domain of T). One can see that (E(X)f,g) is right 
semi-continuous. In applications, the operator T is often a differential operator like 
the Laplacian. Intuitively, this representation generalizes the concept of a singular 
value decomposition of a matrix. 

Following the definition in [7], let any measurable function <j) : (0, 00) — * (0, 00) 
be called an index function. Then 

/>oo 

(/, 5 ) := / cj)(X)dEf. g (X), 11/112 = (/,/), 
Jo 

is a densely defined bilinear form on H with scalar product (•, Let the closure 
of the domain of this bilinear form be denoted H^. Note that then becomes 
a Hilbert space with scalar product (■, and we call the set of all possible H$ a 
variable Hilbert scale. 

In 2 , Begout and Soria introduce a generalisation of the Holder inequality for 
LP spaces. Here, we adapt this theorem to derive a key property of the above 
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variable Hilbcrt scales. In the following, it will be assumed that products of index 
functions are defined point wise; i.e., 

(<^)(A) = <t>(X)9(X). 

Theorem 1 (Generalised Holder inequality). Let <p,ip and 9 be three index func- 
tions and $ and ^ be concave functions (0,oo) — > (0, oo) such that a.e. 

(4) 1<#(#A))*MA)). 
Then 

(5) i<*(ii/iiyn/ii«)*(ii/iiyii/iii) 

for all f e 

Proof. Let the measure v be defined by dv(X) = \\f\\g 2 9(X)dE fJ {X) for / G 
H(<j>+i/>+i)e- By definition, because 9 is positive and the integral of dv equals 1, 
v is a probability measure and, by equation one obtains 

/>oo />oo 

1=/ dv{X)< $ 1 / 2 (0(A))* 1/2 (7/;(A)) ( ij/(A). 



After taking the square of the right-hand side, an application of the Cauchy- 
Schwartz inequality yields 

[>00 [>O0 

1 < / §{(j>(X))dv{X) / V(ip(X))dv(X). 
Jo Jo 

Since both <I> and \& are concave, one can use the (inverse) Jensen inequality to 
obtain 



1 < $ (j 0(A)tfa/(A) 1 * (J V(A)<i^(A) 

The required inequality ([5]) is then obtained by replacing du{X) by its definition. □ 

A related inequality is obtained by choosing 4>(X) — X~ m , ip(X) — X m and 
9(X) = $(A) = \&(A) = 1. In fact, one recovers a special case of a generalised 
Holder inequality which holds for Sobolev spaces (see [TJ p. 50]); namely, 

||/|| 2 < ||/||-m ll/ll™ 

where ||/||_ m := \\fW4, and ||/|| m := 

The interpolation inequality for variable Hilbert scales (see theorem 2.2 in [7]) can 
now be obtained as a direct consequence of the above generalised Holder inequality: 

Corollary 1 (Interpolation inequality 7J). Let (f>, tp and 9 be index functions. 

• // 4> ° ip * s concave, then 

(6) \\f\\%< WfWUo^imlelWfWl), feH^ +l)g , ffLO. 

• If (j) and ip are strictly increasing and ip o <fr x is convex, then 

Proof. Let $(A) = 1/A and *(A) = <f) o if}- 1 (X) . It follows that 

1 



-M- x ^(X)))=l. 



0(A) 

Furthermore, $ is concave, and, as tp° 4'~ 1 is convex and monotonically increasing, 
is concave. The following inequality then follows from Theorem 1 
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which can be rearranged to give the first inequality. The second inequality follows 
from the monotonicity of <\>~ 1 . □ 

The standard Hilbcrt scales inequality is recovered by choosing 0(A) = A a , 4>(X) = 
A b ~ a and tp(X) = A c ~ a with b > a and c > a. If a = (b - a)/(c - a) < 1, then 
<f) o ip^ 1 (X) = X a is concave, and the standard interpolation inequality for Hilbcrt 
scales is obtained 

11/116 <||/||M/llc, 

with \\f\\ b := ||/|| , ||/|| a := ||/|| e and ||/|| c := ||/||^. By choosing 0(A) = 1 in the 
previous corollary, one obtains the version of the interpolation inequality which will 
be used in the sequel. 

Corollary 2. Let cf) and tjj be index functions such that </) o tp^ 1 is concave. Then 

ll/ll^ll/f^dl/liyil/ll 2 )), feHi,. 

A common situation in applications arises when (f>(t) — t and ip(t) = t m for some 
integer m > 1. For m > 1, (f>o -0 _1 (A) = A 1 /" 1 is concave, which, thereby, yields 

n/iii < ii/ir- i/m n/L /ro 

2.1. Dilational Variable Hilbert Scales (DHS). A new special family of 
variable Hilbert scales can be generated from a monotonically increasing index 
function a : [0, 00) — > [l,oo) with a(0) = 1, when the scales are all of the form 
4>(X) — a(sA), with s > 0, and the corresponding norms are defined to be 

a(sX)dE ff (X). 



The Hilbert space with norm || • || s will be denoted by H s . As the index functions 
are obtained from dilations of the original function a(A), we will refer to this family 
of variable Hilbcrt scales as dilational Hilbert scales (DHS). Because the generating 
index function a is monotonically increasing, it follows that, for s < t and / S H t , 
II /lis < || /||t, which implies the existence of a continuous embedding H s H t . 
Since a(0) = 1, it follows that ||/|| = ||/|| and H = H. 

Though ordinary Hilbert scales (OHS) are not directly derivable as a special case 
from DHS, they do form a special subset of DHS as identified in 

Proposition 1. The (ordinary) Hilbert scales, defined by a self-adjoint operator 
T > I, are identical with the dilational variable Hilbert scales generated by a(X) = 
exp(A) and the operator logT. 

Here / is the identity and T > I means that the spectrum is contained in [1, 00). 
This is required so that the logarithm log T defined by 

POO POO 

logT= / log(A)dE(A) = / ^dE(cxp(p)) 
J 1 Jo 

is well defined and positive definite. The proposition follows automatically on 

applying these special choices to a(s^i) = exp(slogA) = A s , where /i — log(A), and 

noting that the spectral projections of log(T) are just E(exp(p)). 

The next theorem identifies an important subclass of DHS which genearte a quite 

special and useful set of interpolation inequality. 

Theorem 2. Let a(X) satisfy the scaling relations 

a(aX) < -a(X), a e (0, 1], 
a 

and, for some c > 0, let the generating function of the DHS be 
a(X) = 1 + c j exp a(s)ds^j dt. 
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Then one obtains, for every a G [0, 1], the interpolation inequality 

WfWlt < ll/llV^di/i^/ii/ii 2 ). 

Proof. This theorem is an immediate consequence of corollary [T] with 0(A) = a(o~X) 
and ip(X) — a(A), once it is established that the conditions of the corollary are 
fulfilled. 

By definition, a(A) is differentiable with derivative 

6(A) = cexp ( J a(s)ds^J . 

Because c > 0, it follows that d(A) > for A > and, hence, that both <f> and tp 
are strictly monotonically increasing. 

As a is an integral with positive integrand, it follows that 0(A) < i/j(X). 

The remaining condition to show is the convexity of ip o _1 . Let 

6(\)=^- 1 (\)) = a(a- 1 (X)/a). 
From the chain rule, one obtains 



0(A) = 



The convexity of 6 is equivalent with the monotonicity of 8, which, as a is mono- 
tonically increasing, is equivalent to the monotonicity of £(A) = #(a(A)). One now 
has 

a(X/a) 



C(A) 

di 

C(A) 



trd(A) 

and so, by the quotient rule of differentiation, 

o(A/ct) 6(A/ct)6(A) 



<t 2 6(A) ct6(A) 2 ' 
Because 6(A) > 0, the required monotonicity reduces to the scaling condition 

a(\/a) > J(X) 



a{\/a) ~ 6(A) 
lows t 
2(A) d 



From the definition of a, it follows that log 6(A) = log(c) + L a(s)ds and, hence, 



6(A) dX ^aW=a(X). 

On replacing A with A/a in the scaling condition for a, the required monotonicity 
and, hence, convexity are established. □ 

The simplest example arises on choosing a(A) = 1 and c = e, which corresponds 
to taking <z(A) = cxp(A). Another example, which is the limiting example for the 
convexity condition, arises when a(A) = (7 — 1)/A and c = 7, which corresponds to 
taking a(A) = 1 + A 7 . Even though this generates a one-parameter family of scales 
like the OHS family , it is a different family because its index function takes the 
form a(sX) = 1 + s 7 A 7 where s is the parameter indexing the Hilbert spaces H s . 

Remark. In the deliberations below, the interpolation inequalities derived above 
proved to be sufficient for the purposes required. Though ordinary interpolation 
inequality can be used to bound the norms ||/||(i-<T)r+<Tt m terms of the norms 
||/||r and ||/||t, for any < r, only the situation with r = is discussed in the 
sequel. Consequently, this represents a matter where the current theory might be 
generalised. In addition, other extensions might include different choices for the 
functions <& and '5 in the generalized Holder inequality. 



DILATIONAL HILBERT SCALES 



AND DECONVOLUTIONAL SHARPENING 7 



3. Solving Ill-posed Problems 

Because of their special structure, variable Hilbert scales interpolatory inequal- 
ities yield a natural framework in which to derive upper bounds for the errors 
associated with approximate regularized solutions of ill-posed operator equations 
Af = g for "non-exact data" g t . For well-posed operator equations, they are not 
required. Assume, for simplicity, that A is injective so that, for every element g in 
the range R(A) of A, there is a unique / such that Af = g. For well-posed A, R(A) 
is the full Hilbert space H, which guarantees that g e € R(A) and, furthermore, 
that the inverse A -1 of A exists and is bounded. Consequently, on assuming that 
1 1 <7 e — g\\ < e, one can immediately derive the upper bound 

\\h-i\\<\\A- l \\\\g e -g\\<\\A- x \\t 

without the need for an interpolatory inequality. 

For ill-posed problems, R{A) is not closed and A -1 , if it exists, is not bounded. 
Consequently, unlike for well-posed problems, there is no guarantee that there ex- 
ists an f e such that the residual r = Af e — g is equal to g e — g. However, from 
the boundedness of A, it follows that \\Af e — g\\ < ||A||||/ e — /||. It therefore fol- 
lows that small residuals is a necessary condition for the errors to be small. In a 
way, this condition corresponds to the consistency condition in the Lax-Richtmyer 
equivalency theorem In fact, for the approximations f e generated by regulari- 
sation methods, the pseudo-data g e := Af e is a well-defined function of f e . One can 
thereby interprete f e as the solution of Af e — g e . Consequently, the consistency 
error corresponds to the difference 

g-g e =Af-g e = g- Af e = -r, 

and, thereby, to the negative of the residual generated by f e . 

However, because the available data are only known approximately as g e , the 
value of the residual Af e — g cannot be determined and utilized to assess the ap- 
propriateness of a given f e . A popular surrogate is the discrepancy Af e — g t . The 
goal is to limit the choice of the f e so as to guarantee that \\Af e — g e \\ < £■ From 
the triangular inequality, it follows then that 

\\Af e -g\\ < \\Af e -g e \\ + \\g e -g\\<2e, 

which is just twice the value obtained for a well-posed problem. 

Interestingly, this bound for the residual, associated with the reconstruction, is 
conceptually similar to the consistency bound used in the theory of the numerical 
solution of well-posed initial value problems and operator equations. In such situa- 
tions, convergence is established by imposing a stability condition on the numerical 
performance of the numerical method and then applying the Lax-Richtmyer equiv- 
alence theorem (see, e.g., For ill-posed problems, the interpolation inequality 
can be used to derive a similar (equivalence) theorem. In fact, the source condition, 
in the form / = Su with \\u\\ < C for some constant C, can be reinterpreted as 
a stability condition by assuming that the reconstructions f t — Su e are such that 
Il u e|| < C f° r some constant C. 

The operator S, which defines the source condition, is often chosen to take the 
form S = (A*A) m , m > 0. On setting the generator T of the Hilbert scales to be 
(A* A)" 1 , the norms in which / are bounded are thereby determined along with 
the associated counterpart of equation (1). However, such source conditions, which 
correspond to taking <p(t) = t a and i/j(t) = t", do not always generate sharp enough 
error bounds. This is resolved by defining S to be some more general function of 
A* A. An example is given in Subsection 3.1, where the differentiation of analytic 
functions is examined. It illustrates how variable Hilbert scales are able to yield 
substantially improved error bounds. 
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If it is assumed that there exists an operator T and index functions </> and tp 
such that R(A) C and R(AT) C H^, one obtains a chain of three embedded 
Hilbert spaces H^, <^-> ff^ <^-> H with g e ff^,. Regularisation is then achieved by 
controlling the size of ||j4/ e ||^. 

Here, it will be assumed that max(||<7||.0, IIAfelU) < C such that \\r\\^ < 2C. One 
then obtains, for e = f e — /, on using the corresponding interpolation inequality 
for variable Hilbert scales (Hegland [7j), the following bound 

||e|| 2 < \\r\\l < \\r\\ 2 cbo^(\\r\\l/\\r\\ 2 ) < ||r|| V ° *~ \W IVt) 

where v — S~ 1 e — u e — u satisfies the bound < \\u\\ + C. 

If it is assumed that (f>(t)/ip(t) is monotonically decreasing such that 

lira </>(*) /V»(t) = 0, 

a convergent rate can be derived. A simple substitution shows that the right-hand 
side of the previous bound is monotonically increasing in ||r||. Consequently, if one 
invokes the discrepancy ansatz (i.e., chooses / e such that ||r|| < 2e and ||r||^ < 2C), 
it follows that 



(8) ||e|| < 2 ev /0o^-i(C 2 e- 2 ) 

which is decreasing monotonically to zero as e — > 0. 

Consider, for example, the case where ip(t) = exp(t) and <f>(t) — t. Because 
= exp(i) is convex and t < exp(t), the conditions given in Hegland [7] 
are satisfied. In this way, the following error bound holds 

(9) ||e|| <2V2(logC/e)e 

for sufficiently small e. 

For DHS, the chain of three Hilbert spaces becomes H t H at H. In 
addition, using the interpolatory inequality for DHS of Theorem^ one now obtains, 
assuming ||e|| < ||r|| CTt , 

llef^llrll^^HrllVao-^llrllf/llrll 2 )). 

For each of the examples considered at the end fo Section it will be shown that 
the assumed inequality |je|| < \\r\\ at is valid. 

3.1. Analytic differentiation exemplification. Ordinary Hilbert scales corre- 
spond to choosing ?/>(A) = A m and <p(X) = A. Because 0o , i/; _1 (A) = A 1 '™ is concave, 
it follows that 

||e|| < 2C 1/m e 1 - 1/m 

which decays much slower to zero than the bound given by @ . Consequently, if it 
is applied to estimating the errors in the computation of the derivative of observed 
data of an analytic function, suboptimal error bounds will result. 

As exemplification, consider the determination of the derivative of a periodic 
function which is the restriction of an analytic function onto the unit circle. In 
practice, this problem can arise when one needs to determine the electric field from 
measurements of an electrostatic potential. We will consider a slightly simplified 
situation below, but the same arguments can be applied to other applications and 
domains. Consider the unit complex sphere B\ C C. Let H 2 (\) be the Hardy space 
of all analytic functions on B\ with the scalar product 

if, 9) = ^ fj W^)g(e 2 ™ 9 )d0, f,g E H 2 (l), 

and, for simplicity, we will also assume that f{e 2 ' K ' lS )d9 = 0, the general case 
can easily be derived from this case. There is a one-to-one correspondence between 
the potentials which are the real parts of the analytic functions / and g. It follows 
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directly that the monomials form an orthonormal system of H 2 (l) and furthermore, 
that the norm in H 2 (\) is 



i/n 2 = E 



k=l 



if f(z) — Yl'kLi a kZ k (recall that we assume that ao — 0). Consider now the self- 
adjoint positive operator T defined by 

oo 

(Tf,g) = E kakh- 

A differentiable function is one where the derivative is in H 2 (l). This requires that 

2 oo 



dg(z) 



dz 



= E fc2 N 2 < 



fc=l 



If we set H = H 2 (l), then the differentiable functions are just the elements of the 
space H,/, defined previously using the operator T above along with = t 2 . 

Here, the source condition can be formulated as a condition on the (unperturbed) 
data g. If it is assumed that there are no singularities of the field in a sphere of 
radius R > 1, then the data g(z) is analytic in the sphere Br and in the Hardy 
space H 2 (R) of analytic functions on Br with a finite squared Hardy norm: 



-, r 2ir oo 

±- / \g(Re 2 ^)\ 2 d6 = Y^R 2k \h\ 2 < 



oo, 



where g{z) = J2T=i^ kzk ■ ^ follows that g e with ip(t) = R 2t . (Actually, it 
is the restriction of g onto the unit sphere B\ that is examined here.). Applying 
the VHS error bound (0 , one obtains the following bound for the differentiation of 
analytic functions 

2e|log(C/e)| 



< 



logi? 



This yields an explicit demonstration that substantially better asymptotic error 
bounds can be derived using VHS than can obtained using OHS. 



3.2. Source conditions. In some cases, the "exact source condition" that gives 
the maximum convergencen rate, i.e., the maximal ip for which Af G may not 
be known. One approach, the Morozov method p~5] then stabilises the norm ||/ e || 
of the solution. A more general method would lead to a bound on ||7"||e for some 
index function 9. In this way, a chain of four Hilbcrt spaces 

is generated. The <j) relates to the error norm which one would like to be small, the 
8 corresponds to the norm which is controlled by the algorithm and the ip relates 
to the smoothness of the solution. In the case discussed in the previous subsection 
for the embedding discussed in the first part of Section [31 one has 9 = ip and, for 
the Morozov method, one has 9 = <fi. For 9 = ip, an algorithm based on ||r||e would 
then yield the bound 

llef^lHlVo^dlrllg/llrf). 
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Consider now the generalised Morozov method where f e is chosen so that ||/ e \\g < 
|| (The Morozov method gets this bound by choosing / c to minimize the 6- 
norm.) Then one obtains 



where the following variant of the Cauchy-Schwartz inequality J\ has been used: 



In the classical Morozov situation, where = <j> and ip = cj) 2 , the following well- 
known y/e convergence estimate holds: 



Combining this with the corresponding interpolation inequality, error rates are 
derived for situations where <f> < 6 < tp. 

3.3. Application - Deconconvolution by sharpening. The theory discussed 
above provides a natural framework in which to analyse the sharpening of broadened 
and possibly overlapping spectroscopic peaks by deconvolution. Because, as shown 
in the introduction, the convolution of two Gaussian peaks is a Gaussian peak, 
an OHS analysis can be applied successfully when the peak, to be measured, and 
its broadening, that occurs as the result of its measurement, are both modelled as 
Gaussians. In the DHS framework, general dilational-parameterized peaks can be 
analysed with similar facility. The VHS framework could be applied, but additional 
assumptions would have to be invoked, like the ones given above in Theorem 2 
for DHS, before error estimates as sharp and useful as those given below could 
be derived. In essence, Theorem 2 generates a framework which allows the VHS 
methodology, as a DHS methodology, to be applied directly to deconvolution by 
sharpening. 

Spectroscopy reveals information about the chemical composition of samples 
and is an important tool in chemistry, physics, biology, astronomy and related 
industrial applications. The data consist of a superposition of "peaks". In the 
case of overlapping peaks, their separation and identification poses a substantial 
challenge. Methods for performing such tasks are discussed in [5]. They have wide 
applicability and can also be used for deblurring in image processing. 

The widening of the peaks in a spectrum results from a "diffusion" of information 
into neighboring frequencies. If this "diffusion" is independent of location, it can 
be modelled as a convolution. In an Z,2(R) Hubert space context, the theoretical 
model takes the form 



where / 6 1>2(R) is the actual spectrum being measured, g € -^(R) is its mea- 
surement and B 1 {x) = r y~ 1 B(x/j) models the broadening that the measurement 
process has performed. The observed spectral data g e are perturbed by "observa- 
tional noise" such that \\g e — g\\ < e. The problem of recovering / from g e is called 
deconvolution. It can be severely ill-posed for smooth B. Because of the importance 
of spectroscopy across a wide range of applications, deconvolution is a prominent 
example of an ill-posed problem. 

Rather than attempting to accurately perform the full deconvolution, it is more 
sensible to "sharpen" the spectrum so that a better identifiability of the locations 



\\r\\ 2 e =\\Aff e -2(Af,Af e ) 9 + \\Af e \ 
<2(g,r) e 



(g,r) e < \\gU IMUw- 
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and number of peaks is achieved compared with that available from a visual inspec- 
tion of the available data g e . As explained in the Introduction, this corresponds to 
finding the solution z of 

B ltj3 * z = g 

for the data g e where the regularity of z is determined by the source condition 

z = Bp*f, /eLa(R). 

When applying the DHS interpolatory inequality in order to obtain bounds on the 
error, it is appropriate at this stage to recast the source condition as a condition 
on the data 

9 = B 1 *f, JeL 2 (R). 
In order to utilise the DHS interpolation inequality, it is necessary to introduce a 
DHS (as introduced in section l2~Tj) which guarantees that g is in H s . One way to 
achieve this is to choose a DHS such that 

(10) \\B y *f\\s = 11/11, feL 2 (R) 

holds for some fixed s. We now proceed to derive a DHS for which this is satisfied. 

For convolutions, it is natural to generate the DHS with T = —d 2 /dx 2 . Using 
the Fourier transform /, one obtains 

XdE ff (X) = -L ^ 2 |/HI 2 du = ^ f™J> (|/H| 2 + |/(-^)| 2 ) 



du 

and dE ff (X) = ± (\f(u)\ 2 + |/(-^)| 2 )- The resulting DHS norm is 
given by 

imi " = h r a{suj2) ( |/h|2 + ,/( ^ )|2 ) duj - 

For the sequel, assumptions must be invoked regarding the Fourier transforms 
B(lu) of the peaks. First, it is assumed that the absolute values of the Fourier 
transforms are symmetric; i.e., that 

\B(u)\ = \B(-u)\. 

A straight forward application of the definition of B 7 (t) and the Fourier transform 
proves that B 7 (lo) = B{^uj). With this and the assumed symmetry, one then 
obtains for the left hand side of condition (1101) 



Since 



2vr_ 

the required condition (jTUJ) holds if s = j 2 and 

1 



du>. 



^ J q (|/MI 2 + |/(-^)| 

slds if s 

(11) o(A) - _ ■ 

\B(VX)\ 2 

We now show that, for three of the most important examples of peaks, this choice 
does indeed lead to a DHS for which Theorem [2] can be applied. For Gaussian 
peaks, one has B(u>) — e~^~ I 2 and, by equation pip , it follows for this case that 
a(A) = exp(A). For this case, the DHS coincides with the OHS. For the second 
example, consider exponential peaks for which B(ui) — 1/(1+oj 2 ). In this situation, 
a(A) = (1 + A) 2 . Going back to the definition of the function a(A) which generates 
a, one finds that a(A) = (1 + A) -1 . It follows that a(a\) < a(X)/a. Consequently, 
such peaks also generate a DHS for which the interpolation inequality holds. For 
the third class, B(uj) — e - '"'. It follows that a(A) = exp(2vA) and hence a(X) — 
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(1 — l/(2\/A))/\/A for which the inequality a(erA) < a(X)/a also holds. Thus, the 
corresponding DHS interpolation theorem holds for the members of this class. 

We have thus seen that, for three of the most important classes of spectra, the 
interpolation theorem holds for the DHS generated. However, the theorem only 
provides a bound for the norm ||r|| scr of the residual r. It is now necessary to use 
it to bound the error norm ||e||. Initially, a condition is imposed on the Fourier 
transform of the peaks which guarantees the bound 

INI < IMU- 

Since r = -B 7l/ 3 * e, it can be shown that, under appropriate conditions, 
INI < \\B y ,p * z\\ s „, z £ L 2 (R). 

As the norm is a continuous function, it suffices to show that this holds for a dense 
subset in L2(R). Furthermore, since z — Bp * /, it is only necessary to prove that 

11-8/3 * /II < II^WIIL- 
From the Fourier transform, it is clear that this inequality holds when 

\B(^)?< ^ (7 " } ' 2 , «6E 
" |B( 7 ^)| 2 

Some rescalings reveal that there exists a a < 1 for which this holds if and only if, 
for each scalar < r < 1, there exists a < p < 1 such that 

\B(tu)B(pw)\ < \B(w)\, uj G R, 

or, in terms of the functions a defining the DHS, the condition becomes: for every 
r, there has to exist a p such that 

(12) o(A) < a(rX)a(p\), A > 0. 

In the first example above, since a(A) = e A , equality in (|12p is guaranteed if t + p = 
1. In the second example, since a(A) = (1 + A) 2 and, because (1 + tw)(1 + pu) > 
1 + (t + p)uj, the inequality holds whenever r + p = 1. In the last example, since 
a(A) = exp(2v / A), equality holds if yfr + y/p =1. It follows that, for the peaks 
considered, the choice of DHS suggested above does allow for the application of the 
interpolation inequality. In addition, the interpolation inequality provides a bound 
for the error associated with the reconstruction procedure. Of course, the peak 
sharpening procedure does need to satisfy the conditions discussed in the previous 
sections. 

For the Gaussian, exponential and rational peaks, the conditions required for 
the application of the DHS interpolation inequality are satisfied. For the resulting 
bounds on the norm of e = z e — z, one obtains 

(i) the usual OHS error bound for the case of Gaussian peaks where a(A) = e A : 

Ml < IM^Mi 1 -; 

(ii) a convex combination for the case of exponential peaks where a(A) = (1 + 
A) 2 : 

||e||<a||r|| s + (l-<7)||r||; 

and 

(iii) for the rational peak where a(A) = cxp(2VA): 

Hell < IMI^IIrll 1 -^. 
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In summary, one obtains from the interpolation inequality three types of "convex 
combinations" for the errors, with the a determined by the reconstruction method. 
It follows from what was said above that a satisfies a + f3/j = 1 for the first two 
shapes and ^fa + \[Wfl = 1 f° r the third. By tuning (3, one thereby tunes the 
parameter a. Ideally, one would like to choose er = to get the smallest possible 
error. This, however, corresponds to the case (3 = 7; i.e., to the case where no 
sharpening is done. Conversely, (3 — corresponds to "full sharpening" but in 
such situations, the errors will be large. In the practical application of the above 
bounds, it is necessary to carefully consider how the trade-off is performed between 
the amount of sharpening and error in the reconstruction of the sharpened peaks. 
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